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Abstract. - Critical curves "critical current - external magnetic field" of long Josephson junctions 
with inhomogeneity and variable width are studied. We demonstrate the existence of the regions 
of magnetic field where some fluxon states are stable only, if the external current through the 
junction is different from zero. Position and size of such regions depend on length of the junction, 
its geometry, parameters of inhomogeneity and form of the junction. The noncentral (left and 
right) pure ffuxon states are appeared in the inhomogeneous Josephson junction with increase 
in the junction length. We demonstrate new bifurcation points with change in width of the 
inhomogeneity and amplitude of the Josephson current through the inhomogeneity. 
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Introduction. — Starting from the classical paper 
of Owen and Scalapino, [1] the long Josephson junctions 
(•LJJ) are attracting attention both experimentalists and 
theoreticians. They are good candidates for a wide range 
of applications such as superconducting quantum interfer- 
ence devices, Josephson voltage standards, logic elements, 
and Josephson flux-flow oscillators because of their non- 
Ifnear behavior and quantum effects. [2-4] From the other 
side, the nonlinearity makes the LJJ a very complex sys- 
liem, and some aspects of its properties are not investi- 
gated up to now. Particularly, the experimental methods 
1^0 determine the vortex structure under the critical curve 
are not developed well yet. But there are different kind 
of stable states in the LJJ, which bifurcation curves are 
close to the critical curve of the junction and manifest 
themselves in the experiments. [2, 5] Particularly, some 
"metastable" states (experimental points under overlap 
curve) were observed in Ref. [6]. Static bound states of 
fluxons in LJJ with artificial inhomogeneities were exper- 
imentally found in Ref. [5]. Experimental and theoretical 
investigation of such states is an interesting and actual 
problem. 

An important step of the sample characterization in 
the experiment with LJJs is the measurement of its static 
properties, i.e. the dependence of the critical current Ic 
on external magnetic field H . This dependence provides 



a technique to evaluate the several important parameters 
of LJJ from the experiment, particularly, a critical cur- 
rent density, a magnetic flux penetration field, an effec- 
tive magnetic thickness. [7] Since the solution of the cor- 
responding boundary value problem (BVP) cannot be per- 
formed analytically in general case, the most straightfor- 
ward way to study the dependence Ic{H) is its numerical 
simulation. 

In our previous study we investigated the stability of 
the vortices in LJJ with inhomogeneities (LJJI) [8, 9] and 
vortex structure in exponentially shaped Josephson junc- 
tions (ESJJ). [10-12] The obtained bifurcation curves of 
the mixed fluxon- antifluxon states in LJJI [9] and pure 
fluxon states in ESJJ demonstrate the intervals of mag- 
netic field, where these states are stable only, if the current 
through the junction is not equal to zero. These peculiari- 
ties were not mention at that time. We call these states as 
the CbC-states (created-by-current), and the correspond- 
ing regions of the external magnetic field as CbC-regions. 

In this Letter we study the CbC-states of mixed fluxon 
states $~^<i>^, and pure fluxon states We 

demonstrate new bifurcation points with change in width 
of the inhomogeneity and amplitude of the Josephson cur- 
rent through the inhomogeneity. The influence of the 
model parameters on CbC regions in the LJJI and ESJJ 
is investigated. 
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Method of calculation. First we start with LJJI 
case. We consider the Josephson junctions in the overlap 
geometry. [2] In order to study a hnear stability of the 
static distribution ip{x) of magnetic flux (j}{t,x)^ we write 
the perturbed expression for x) as x) — ip{x) + 
e^{t)'ip{x), where e is a parameter [13, 14]. Then in a 
first approximation with respect to e, we get the following 
boundary value problem (BVP) for (p{x) 

-(Pxx+jc{x)smip + -f = 0, ip{-l) = ip{l) = he (1) 

and the corresponding Sturm-Liouville problem (SLP) 

-ipxx + cr{x)ipx + q{x)tp = Xi^ i'xi-l) = Ml) ^ (2) 

Here all the quantities are in diniensionless form [3] , is 
the external magnetic field, 7 is the external current, and 
A is an eigenvalue of the SLP, 21 is a length of the junction. 
The function jc (x) represents the Josephson current am- 
plitude, so I jc(2;)| < 1- Because of the nonlinearity the 
BVP, eq.( 1) has more then one solution for given param- 
eters (counted set in case /le = 0, 7 = 0, Z — > 00) [15]. The 
potential q{x) = j oix) cos ip{x) is defined by the concrete 
static solution ip{x), which is asymptotically stable with 
respect to small perturbations, if lim \£,{t)\ ^ 0. As the 

t — >oo 

junction's length / < 00 in case under consideration and 
the potential \q{x)\ < 1, the SLP eq.( 2) is regular, so 
there exists a discrete low-bounded set of real eigenvalues 
An, n = 0, 1, 2, . . . and the minimal one Aq > —1. [16] The 
corresponding eigenfunctions satisfy the norm condition 

/ '4P'{x)dx — 1 = 0. In case Aq < the static solution 

^p{x) is unstable, i.e. | 00, when t — > 00. Solutions 

of the BVP eq.( 1), eq.( 2) describe the different fluxon 
states in LJJ. Numerical simulation simplifies the study 
and makes it possible to estimate the range of variation of 
the parameters in which one can expect stability or insta- 
bility of the magnetic flux distributions in the Josephson 
junction. 

The models of LJJ depend on m parameters, such as 
junction's length 21, external magnetic field he, external 
current 7 and so on. Let us denote the vector of all the 
parameters by p. In all cases we presume that the possible 
solutions of non-linear BVP eq.( 1) continuously depend 
on the set p, i.e. Lp = ip{x,p). It follows that the potential 
of SLP eq.( 2) generated by the solution ip{x,p), depends 
on the model parameters p as well, i.e., q = q{x,p). Fi- 
nally, the corresponding eigenvalues and eigenfunctions of 
eq.( 2) are continuous functions of 77, i.e. A„ = A„(p), 

i^n = 1pn{x,p). 

Hence, the static solution tp{x,p) is asymptotically (ex- 
ponentially) stable with respect to small space-time per- 
turbations in some bounded subset of the parameter re- 
gion, if the minimal eigenvalue satisfies Xo{p) > 0. [14] 
Every point on the hypersurface Aq (p) = is the bifurca- 
tion point for the solution under consideration. The values 
of the parameters which satisfy Xo{p) = 0, are called the 



bifurcation (critical) values for this solution. The cross 
section of this hypersurface by a hyperplane, which cor- 
respond to fixed values of m — 2 parameters, determine 
the bifurcation curve for the rest two parameters. From 
the experimental point of view, the most important are 
the bifurcation curves of kind "critical current-magnetic 
field" 

Xn{j,he) = 0. (3) 

The numerical algorithm for the determination of the bi- 
furcation curves is proposed in [17] and was successfully 
applied to various physical problems [9] 

Bifurcation curves in LJJ with inhomogeneity. 

— The static configurations of finite-length junction in 
the presence of an external magnetic field were studied by 
many authors. [1,8,9,14,18-21]. Particularly, S.Pagano et 
al. [18] for arbitrary junction lengths demonstrated the bi- 
furcation curves of the fluxon states in the three different 
regimes typically observed: small, intermediate, and long 
junction. We consider here the LJJ which barrier layer 
contains one resistive rectangular inhomogeneity, charac- 
terized by its width A, the position of the center of the 
inhomogeneity C and the portion of Josephson current kj 
through it. The existence of the inhomogeneity leads to 
the local change of the Josephson current through the 
junction, which can be modeled by jci^) = 1 + fcj in- 
side of the inhomogeneity and jc{x) = 1 outside of it. 
At kj > the value of the current through the inhomo- 
geneity exceeds the value of the current in other parts of 
the junction and such inhomogeneity are considered as a 
shunt. [14] At kj £ [—1, 0) the value of the current through 
the inhomogeneity is less than in other parts of the junc- 
tion and such inhomogeneity represents a microresistor. 
Changing the thickness of the barrier layer inside of the 
inhomogeneity, we can model the transformation from the 
"shunt" to the "microresistor" inhomogeneity. The value 
kj = corresponds to a homogeneous junction. Results 
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Fig. 1: Critical curve of the LJJ with length 21 — 7 and width 
of inhomogeneity A = 0.7. The insert shows the bifurcation 
curves in the marked by circle region. 

of numerical solution of non-linear eigenvalue problem [9] 
eq.( 1), eq.( 2) are presented in fig. 1, where the curve 
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Fig. 2: (a) The distribution of the magnetic field tpx{x) along 
the junction for bifurcation states M, and in LJJ 

with length 21 = 7, width of inhomogeneity A = 0.7 at fee = 
1.9; (b) The distribution of the magnetic field for LJJ with 
21 = 12 for $\ ^] and $^ at he = 1.5 and 7 = 0. 



"critical current - external magnetic field" of the junction 
with 21 = 7 and the inhomogeneity with width A = 0.7 
in the center of the junction is shown. This curve is "con- 
structed" as envelope of bifurcation curves eq.( 3), which 
belong to different possible bound states [9]. The inho- 
mogeneity leads to the appearance on the critical curves 
the fragments of bifurcation curves of the stable mixed 
states like $"<()-" and n = 1,2,..., which are 

not stable in homogeneous case. In one's turn this leads 
to the non-monotonic decrease of maximums of critical 
curve, when the magnetic field increases. As we can see in 
fig. 1 , in the investigated region of magnetic field the criti- 
cal curve is an envelope of the critical curves for the Meiss- 
ner state M, mixed distributions <I>^$~\ $-2(j)2^ $-3^3 
and pure fluxon states <i>^, <i>^, <I>'^. All "even" pure 
states like 3>^, etc., are unstable, so their bifurcation 
curves eq.( 3) do not share for critical curve of junction as 
whole. 

The insert to the fig. 1 shows the details of formation of 
the critical curve in the marked region. In all our figures 
in this paper the solid lines correspond to the positive 
direction of current, the dashed lines correspond to the 
negative direction of the current. The critical curve of the 
junction consists of pieces of the bifurcation curves for <I>^, 
$^<I>^^ and states in this region. 



Table 1: Physical parameters for bifurcation states M, 
$-^$1 and $i<l>-i 



Type 


7cr 


N 






$1 


0.228 


1.317 


1.69 


1.438 


M 


0.065 


0.073 


0.847 


0.048 




0.034 


1.654 


2.28 


1.58 




0.017 


1.43 


2.19 


1.329 




-0.017 


2.57 


2.19 


2.671 




-0.034 


2.346 


2.28 


2.42 


M 


-0.065 


-0.073 


0.847 


-0.048 


$1 


-0.228 


0.683 


1.69 


0.562 



The bifurcation curves for mixed $1$"^ and $^1$! 
states demonstrate an interesting peculiarity. As we can 
see, in the intervals of magnetic field {hi, ^12) and (/13, /14) 
these CbC-states are stable only if the current through 
the junction is not equal to zero. A value of the cur- 
rent, which is needed for the creation of stable mixed 
distributions $i<I>^i and <I>^i$i in CbC-regions, depends 
on the value of the external magnetic field. The table 
presents the values of the physical parameters for bifur- 
cation states and at /i = 1.9 in compar- 
ison with AI and $^ . Here the average value of the 

magnetic flux (number of fluxons) is defined as a func- 
i 

tional N[(p] ~ ^ J (p{x) dx and the full magnetic flux 

-I 

through the junction is A(p = (p{l) — if{—l). [9] For non- 
bifurcation pure fluxon states at 7 = the number of 
vortices equals to the integer number: A^[<I>='="] = n, while 
for non-bifurcation mixed fluxon states at 7 = we obtain 
7V[$"$-"] + iV[$-"$"] = 2n. [9] For non-bifurcation sta- 
ble Meissner solution iV[M] = 0. As we can see from the 
Table, the found bifurcation values satisfy the following 
relations. For $" states the half of sum of N for posi- 
tive and negative currents equal to n. For CbC states we 
should take the sum of iV[$"$-"] and Af[$-"$"] for op- 
posite current directions to have integer number 2n. The 
magnetic flux in the center of junction Lp{Q)/2'r: fulfil the 
similar relationships. 

The distribution of the internal magnetic field ^x{x) 
along the junction for bifurcation states M, 
and <I>~i$^ at the different values of 7cr, presented in the 
Table is shown in fig. 2(a) in case of LJJ with 21 = 7 
and A ~ 0.7 at /i = 1.9. The inhomogeneity in LJJ at- 
tracts the bound states and it leads to the deformation of 
the distributions of the magnetic field for these states. In 
fig. 3 we present the bifurcation curves for vortices in LJJI 
with length 21 = 12 and the same width of inhomogeneity 
A = 0.7 in the center of junction and the zero amplitude 
of Josephson current through inhomogeneity jc = 0. It 
shows clear the tendency of changing of the different bi- 
furcation curves with increase in magnetic field. The main 
features of this critical curve coincide with the results for 
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Fig. 3: Bifurcation curves for LJJI at 21 — 12 with the inho- 
mogeneity in the center with A = 0.7 

LJJ with 21 = 7. Particularly, we observe here the mixed 
states as well, and as a result, the non-monotonic decrease 
of maximums of the junction's critical curve with magnetic 
field. 

But in contrast to the case 21 ~ 7, the new "noncentral" 
fluxon states are appeared here in addition to the pinned 
to the inhomogeneity "central" ones. We found that in 
case under consideration the noncentral bound states , 
^l, and <I>~^ are stable only. The distribution of the in- 
ternal magnetic field along the junction for <i>;^ and <I>^ 
is shown in fig. 2(b). The new moment in case of 21 = 12 
is a fact, that the states and <I>^ are not stable without 
electric current, at any value of external magnetic field. 
We have not found the stable $1$"^ and states 
in this case also, so we consider that they are replaced 
with and <1>~^ when the length of LJJ is increased. Or 
we can say, that for LJJ with 21 = 12 the width's value 
A = 0.7 is small enough to stabihze the and $"^$1 

states. We can see below in fig. 5(a) the transformation of 
the state into and with the increase in 

A. The $f and demonstrate the CbC- regions which 
look like CbC-regions for mixed states <I>i$^i and $-i<l)i 
for small junction length. Fig. 3 allows us to represent 
from the unified point of view the influence of the inho- 
mogeneity on the bifurcation curves of the fluxon states in 
LJJ. It leads to: a)mixed fluxon- ant ifluxon states 
b)non-monotonic decrease of maximums of the junction's 
critical curve with magnetic field; c) right and left fluxon 
states; d)diffcrent position of the CbC-states of the differ- 
ent fluxon states. 

The influence of the inhomogeneity vifidth. — 

The influence of the inhomogeneity width on the CbC- 
regions for mixed state $1$"^ in the case of 2/ = 7 is 
demonstrated by fig. 4. It shows the dependence of the 
minimal eigenvalue Aq on the external current 7 for <I>^ 
and <I>i(I>~^ for the LJJ with two values of width of inho- 
mogeneity A = 0.5 and A = 1 at /le = 2.2. The zeroes of 
Aq determine the critical currents for corresponding states. 
The two values of current ^cr determine the down and up- 



per critical currents for the CbC states at fixed value of 
the external magnetic field he = 2.2. We can see that with 
increase in A the interval of the CbC-region for $^$-1 
is increased. The influence of the width of the inho- 




External current 

Fig. 4: Dependence of the minimal eigenvalue Xmin on the 
external current 7 for and $^<I>~^ at two different values of 
the width of inhomogeneity A = 0.1 and A = 0.5. 

mogeneity A on the bifurcation curve for is demon- 
strated in fig. 5(a). We show here the bifurcation curves 
at A = 0.3; 0.5; 1. This figure demonstrates the transfor- 
mation — > (4>^$~^, $~^$^) with increase in A. There 
is a bifurcation region in the interval 0.3 < A < 0.5, where 
the state $^ is getting unstable, but <l)i(i)-i and $~^4>^ 
are stable. The details of such transformation are shown 
in fig. 5(b). As we can see, the state $^ is stable yet at 
A = 0.45, but the corresponding region of magnetic field 
decreases with A. This region disappears completely in 
the interval 0.45 < A < 0.5. 

Actually, fig. 5 shows in the current versus mag- 
netic field plane typical features of the catastrophe the- 
ory. [22] This correspondence is stressed by fig. 5(c), where 
the dependence ^(7) of the total energy on 7 is shown for 
stable distributions <I>^, <I>~i$i and related to them un- 
stable ones at he = 1.85 (see fig. 5b). The cusp in the Bi 
point is a bifurcation point for with change in the ex- 
ternal current. Correspondingly, the points B2 and B3 are 
the cusps of the bifurcation points for mixed distribution 

Note, that the interval along the 7-axis which cor- 
responds to the points B2 and B3 has not the point 7 = 0, 
i.e. the current in B2 point is a "creation current", and in 
the point B3 - "annihilation current" of the distribution 

The states the two horns (in the interval of mag- 
netic field hi < h < h2 in the insert to fig. 1) reflect the 
fact that these states can not be stable without current 
through the junction. The detailed study of the observed 
features in the correspondence of the catastrophe theory 
will be done somewhere else. 

Shifting of the inhomogeneity. Fig. 6 shows the 
bifurcation curves for <I>^ in LJJI with 2/ = 7, A = 0.7 
and kj = —1 at different position of the inhomogeneity 
C = 0; 0.5; 1; 1.5; 2 and 2.5. We found that when the inho- 
mogeneity is shifted from the center of the junction, the 
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Fig. 5: (a) The transformation of the critical curve for "1>^ 
with the value of inhomogeneity's width; (b) The details of 
the transformation of $^ at A = 0.45; (c) The dependence of 
the total energy F on external current 7 for distributions 
$-^'1'^ at he = 1.85 
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Fig. 6: Bifurcation curves for $ at different position of the 
inhomogeneity. 



pure fluxon states have the CbC-regions as well. The case 
= correspond to the position of the inhomogeneity in 
the center of the junction. 

The influence of the current value through the 
inhomogeneity. — The decrease in the value of the pa- 
rameter kj, which characterizes the amplitude of Joseph- 
son current through the inhomogeneity, leads to the same 
transformation <3>^ ($^$~^, $~'^$^) as we have ob- 
served with increase in A . We demonstrate it for LJJI 
with 21 = 7 in fig. 7(a), where the transformation of the 
dependence N{he) for <i>^ with the value of the amplitude 
of Josephson current through the inhomogeneity is shown. 
As we can see from this figure, a bifurcation of the states 
exists in the interval — 1 < fcj < —0.735. We found that 
a new bifurcation point exists at kj « —0.77. The de- 
pendence N{he) in LJJ with 21 7, A ^ 0.7, C = and 
kj = —1 is shown in fig. 7(b). It demonstrates clear the 
relations for N for the pure and mixed bifurcation states 
we mentioned above. Particularly, for pure fluxon state <f>'^ 
sum of the N{he) for positive and negative current equals 
to 3, and for states $-2(^2 Ar[$2<j)-2] _|_ 
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Fig. 7: a) Transformation of the dependence N(he) for with 
the amplitude of Josephson current through the inhomogene- 
ity; (b) The dependence N{he) in LJJ with 2Z = 7, A = 0.7, 
= and kj = —1. 

Exponentially shaped LJJ. — As a second example, 
we consider an exponentially shaped LJJ [10]- [12]. In 
Ref [12] we demonstrated that the ESJJ is equivalent to 
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Fig. 8: (a) Bifurcation curves for ESJJ at 2Z = 7 and a = 0.07; 
(b) Bifurcation curves for "I"^ at different values of a. 



the Josephson junction with distributed inhomogeneity. 
So we may expect the CbC regions in the critical curves 
of the ESJJ as well. The corresponding BVP for the ESJJ 
has the form: 

-(fixx + sin ip + a{ip^- he) -J = O^ipxi^l) = K^tpxil) = K 

(4) 

where < cr <i; 1 is the form parameter. 

Results of numerical solution of this BVP combined 
with the corresponding SLP [10] are presented in fig. 8(a), 
where the bifurcation curves for Meissner and four first 
fluxon states are shown for / = 7 and a = 0.07. As 
we can see, the bifurcation curves of pure fluxon states 
in ESJJ demonstrate the CbC-regions. The CbC-states 
are realized at positive current in low magnetic fields and 
negative current in high magnetic fields. The exponential 
shape of the junction leads to the appearance of the dis- 
tributed "geometrical" current a {(px — he), which modifies 
the bifurcation curves of the states in compare with the 
rectangular homogeneous junction. [12] Bifurcation curves 
for <I>^ at Z = 7 and different values of a are shown in 
fig. 8(b). The increase in a leads to the increase of the 
CbC-region. With increase in cr a situation might be real- 
ized, when the dependence of the upper critical current on 
the external magnetic field in CbC regions is getting non- 
monotonic. We found that such non-monotonic behavior of 
the bifurcation curve in low magnetic fields appears in 
ESJJ with increase in junction's length as well. 

Ill summary, we performed the numerical simulation of 
the critical curves of long Josephson junctions with inho- 



mogeneity and variable width. In both cases we demon- 
strate the CbC-regions of magnetic field, where fluxon 
states are stable only if the current through the junction 
is different from zero. We showed that the position and 
size of these regions depend on length of the junction, its 
geometry, parameters of inhomogeneity and form of the 
junction. We consider that development of the experi- 
mental methods for detection of CbC-states will open a 
perspective for their applications. 
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